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Abstract 

This paper presents an algorithm for computing Groebner bases based upon labeled poly- 
nomials and ideas from the algorithm F5. The main highlights of this algorithm compared 
with analogues are simplicity both of the algorithm and of the its correctness proof achieved 
without loss of the efficiency. This leads to simple implementation which performance is in 
par with more complex analogues 
00 . 

£N| ■ Consider polynomial ring P — k[x\, . . . , x n ] over field k. Also assume that monoid of its 

monomials T has a monomial order -<. A problem asking for a Grobner basis can be stated for any 
ideal(/i, ...,/;) in this ring. One of the approaches to the problem is using iterative method which 
computes every step a basis for ideal (/i, ...,/,;), i — 2 . . .1 based on the already computed for 
(/i, . . . , fi—i) basis Ri-i and polynomial /j. The algorithm described in this paper is designed to 
perform one step of such computation. So, the algorithm's input data consist of a some polynomial 
/ and a polynomial set referred as {g\, . . . ,g m } which is Grobner basis of ideal Iq = {g\, . . . ,g m )- 
After finishing the algorithm should give the resulting polynomial set R being a Grobner basis 
of ideal I = (gi, . . . , g m , f). The special cases / = =>• / = Iq and 3i gi E k I = P are not 
interesting from the computational point of view, so the further chapters assume that / ^ 0, \/i gi £ 
k. The homogeneity of input polynomials is not required unlike the F5 algorithm described in [3]. 

o i 

^ ■ Definitions 

VO ■ 

• . Consider the setTo = T U {0} - the monomial monoid extended by zero. The order -< can be 

extended to To as -<o with definition \/t E T t >-q which keeps the well-orderness property. The 
^v^j . notion of division also can be extended to To: ii|^2 — =^3 £1^3 = £2- For polynomial p£P,j)^0 

the highest by -< monom and coefficient are written as HM(p) 6 T and HC(p) £ k. For zero we 
define: HM(0) = £ T , HC(0) = G k. The least common multiple of h,t 2 6 T is written as 
LCM(£i , £2) G T. In the following all definitions are given for fixed Iq and /: 



Definition 1. The labeled polynomial is a pair h = (o~,p) G Tq x P, that satisfies the correctness 



property: 3u G P HM(tt) = a,uf = p (mod Iq). Some terminology is extended to labeled polyno- 
mials. The highest monomial is HM(h) = f HM(p) and coefficient is HC(/i) = HC(p). Additionally 
the signature is defined S(h) = a and a notation is introduced for the polynomial - second ele- 
ment of pair: poly(/i) = p. The set of all labeled polynomials is written as H C To x P. The 
trivial examples of labeled polynomials are (1, /) and (0, g) for g G Io- Another labeled polynomial 
example is (HM(g),0) for g G Iq. It satisfies correctness property because we can take u equal to 
9- 

Lemma 2. The product of h G H,t G T defined as th = (ta,tp) G H, is correct. 
The correctness property is checked by directly finding u for th. 
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Definition 3. If h[, h 2 G H,t £ T satisfy S(h[) y 5(tt 2 ),HM(/i'i) = HM(tft 2 ) ^ 0, then exists 
a signature-safe reduction h[ by h 2 , resulting in labeled polynomial hi £ H, equal to: 

hi = {S(h[),j>olj(h[) + Ktpoly(h 2 )) , 

where the K £ k is selected in a way to perform cancellation of high coefficients, so we have 
HM(/ii) -<o HM(ft'i). Such reduction is equivalent to plain reduction with high term cancellation 
extended with requirement for reductor's signature being smaller than the signature of labeled 
polynomial being reduced. Like in previous case the correctness check is performed directly. 

Let's introduce a partial order <h on H: 

hi = (cri,pi) < H h 2 = {cr 2 ,p 2 ) U HM(pi)cr 2 HM(p 2 )(7i. 

The elements with zero signature or zero high monomial are extremums: 

y<7i,a 2 ,pi,p 2 (0,pi) ftn (a 2 ,p 2 ), (a i,0) ^ H (o- 2 ,p 2 )- 

Lemma 4. Let h\, hi £ H, t £ T. Then hi >h h 2 <^> hi >h th 2 . 

Deduced from the fact that multiplying one of the compared labeled polynomials by t leads to 
multiplying by t both sides in the definition of>n- 

Lemma 5. Let hi,h 2 £ H, BM(hi)\ HM(h 2 ), KM(h 2 ) ^ 0. Then signature- safe reduction h 2 by 
hi is possible ifjhi >h h 2 . 

Deduced from the fact that claims of both sides are equivalent to S(h 2 ) >~o S(hi) ■ 

Lemma 6. Let hi £ H be a result of signature-safe reduction of h[ by some other polynomial. 
Then hi <h h[. 

Deduced from equality S(hi) — S(h[) and decreasing HM during reduction: HM(fti) -<o 
HM(h'i). 

Lemma 7. Let hi <h h 2 be labeled polynomials. Then for W13 € H \ {(0, 0)} at least one of the 
following two inequalities holds: hi <h h 3 or h 3 <h h 2 . 

The lemma clause gives inequality 

HM(fci) S(h 2 ) ^0 HM(h 2 ) S(hi) (1) 

which shows HM(h 2 ) ^ Q,S{hi) ^ 0. Therefore for the special case HM(h 3 ) = we get /13 <h 
h 2 &i\& for the case SQ13) = we get hi <h ^-3- For remaining generic non-zero case the inequality 
((T|) can be multiplied by non-zero monomial HM(h 3 ) S(h 3 ) : 

RM(h 3 ) S{h 3 ) RM(hi) S(h 2 ) ^0 HM(^ 3 ) S(h 3 ) EM(h 2 ) S{hi). (2) 

So, the element HM(h 3 ) 2 S(h 2 ) S(hi) £ To need to be >~o than left side or -<o than right side of 
inequality @, and gives after cancellation one of the inequalities from the lemma statement. 

Algorithm 

Input: polynomial set {gi, . . . , g m } being a Grobner basis; polynomial /. 

Variables: R and B - subsets of H; (o-,p') £ H - current step's labeled polynomial before reduc- 
tion; (a, p) - the same after reduction; r, b - elements of R and B 

Result: Grobner basis of ideal/ = {gi, . . . , g m , f) 
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SimpleSignatureGroebner({(7i, . . . ,g m } , f) 

1. R <- {(HM(ft), 0) , (HMfe), 0) , . . . , (HM( 5m ), 0) , (0, 9l ), (0, <? 2 ), . . . , (0, g m )} 

2. B^{} 

3. (tr,j/)<-(l,/) 

4. do forever: 

(a) p <— ReduceCheckingSignatures(<7,p', B) 

(b) B^-BU{(a,p)} 

(c) if p ^ 0: 

i. for {reR\r < H (<7,p) , HM(r) ^ 0}: 
A. ^Buf ^WyW' r} 

ii. for {r G B | r >h 

A. ^EU{ ^M (g>p)} 

(d) B <- B \ {6 6 B | 3r e Br < H 6 A 5(r)| 5(6)} 

(e) if B 7^ 0: (cr,p') ^— element of B with -<-minimal signature 

(f) else: break 

5. return {poly(r) | r € R} 

ReduceCheckingSignatures(cr, p, R) 

1. do while 3r e Rr > H (cr,p) A HM(r)| HM(p): 

(a) p 4— signature-safe reduce p by >H-maximal element r from the set in cycle clause 

2. return p 

Lemma 8. All pairs from To x P appeared in the algorithm are labeled polynomials fromH\{(0,0)}. 

The elements created before entering main cycle are labeled polynomials mentioned above as 
examples. All other labeled polynomials in the algorithm are created either with multiplication by 
t € T or with signature-safe reduction, so they satisfy the correctness property and belongs to H. 

The clauses of cycles extending B enforces the absence in B elements with zero signature or 
zero highest monomial. So, a never can be and the only R elements with zero signatures are 
(0, <7i), (0, g m ). Any labeled polynomial added to R can have zero highest monomial but R does 
not contain zero polynomial with zero signature. 



Algorithm termination 

Lemma 9. At the any moment during the algorithm execution any labeled polynomial from B can 
be signature-safe reduced by some element of R. 

Labeled polynomials are added to B in a way ensuring existence at least one possible signature- 
safe reductor. The pair (tr, p) G R is such reductor for polynomials added in first for cycle, and 
r £ R - for the polynomials added in the second cycle. 

Lemma 10. Before reduction of polynomial p' - at the ste f[4a\ of any algorithm iteration - the 
signatures of elements{r <E R\r <h (a,p')} does not divide a. 

This holds at the first algorithm iteration because a = 1 and R does not contain elements with 
signatures dividing 1. This holds during next iterations because the existence such elements in R 
would lead to removal (cr,p') from B during previous iterations at the step I4dl 
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Lemma 11. After reduction of p' to p - at the step 4b of any algorithm iteration - the highest 



monomials of elements {r G R \ r >h (cr,p)} does not divide HM(p). 



The cycle in the ReduceCheckingSignatures((7,p, R) stops only when it achieves p for which 



there is no such elements in R . 

Lemma 12. After reduction of p' to p - at the step^4j^ of any algorithm iteration - no one from 
R elements has simultaneously a highest monomial dividing HM(p) and a signature dividing a. 

The lemma IHl ensures that p'is reduced at least once, so (cr,p') >h (o~,p)- Now, by lemmg[7]for 
Vr G R we have r >h (c, p) or r <h {o~,p')- These inequalities allow to apply either lemma [TU] or 
lemma 1111 



Theorem 13. The algorithn SimpleSignatureGroebner({gi, . . . , g m } , /) terminates 



To prove the termination we need to show that all do cycles stops after finite number of 



executions. In the cycle inside ReduceCheckingSignatures(er, p, R) with non-zero p during every 



iteration we have HM(p) decrease according to -<o, which is possible only finite number of times. 
When p becomes zero it stops immediately because of <H-minimality of (<r, 0). 

The set R C To x P is extended every step of the main algorithm cycle. It can be splitted 
to i?« U i? * U where R* C T x {0} , i? * C {0} x P \ {0} , C T x P \ {0}. i? * does 
never extend because a 0. For sets i?*o and i?** we apply a method based on idea of monoid 
ideal introduced in [8] as "monoideal". Consider the following two sets which are monoideals: 
L*o = {{o- 1 (cr,0) G i?* }) C T and L** = ({(tr,t) | 3(a,p) e = HM(p)}) C T x T. Lemma 
1121 shows that elements being added to R expand either L*o or L** in every cycle iteration. The 
monoids T and T x T are isomorphic to N n and N 2ra , so the Dickson's lemma can be applied to 
their monoideals. It states exactly the needed fact - only finite number of expansions is possible 
for such monoideals. 



Correctness of output 

Definition 14. S-representation of h G H over set {r^} C H is an expression poly(/i) = 
]T\ Kjtjpolyiry), Kj G k,tj G T,ij G N, such that Vj HM(/i) ^ HM(t 3 r h ),S(h) ^ S(t 3 r h ). 

Lemma 15. Let poly(/i) = ^ - Kjtj poly(r^) be S-representation ofh. Then at least one j satisfies 
HM{7>.) IIMi/,r, i. 

To get a j satisfying the equality we can take a value which gives the ^-maximum of HM(tj-rj.). 
The next definition extends the notation of S-basis from [T] : 

Definition 16. We call a labeled polynomial set R C H S-basis (correspondingly S a -basis), if all 
elements of H (correspondingly {h G H \ S{h) <q a}) have S-representation over R. 

Lemma 17. Let a 0, R = {rj} be S a -basis and h\,h-z G H,S(hi) = a be labeled polynomials, 
that can't be signature-safe reduced by R elements. Then HM(/ii) = HM(ft,2) and h\ has an S- 
representation over R U {hi}- 

We have from the definition of H that 3ui G P HM(iti) = a, uif = po\y(hi) (mod Iq), i = 1, 2. 
It means that there exists a linear combination of poly(/ii) having signature o~. This can be 
written as: 

3K ek,v e P EM(v) = a' -< a,vf = poly(/ii) - Kpo\y(h 2 ) (mod I ), 

or in the terminology of labeled polynomials: (cr',p') — (a', poly(/ii) — K poly(/i2)) G H. From 
the definition of So-basis and the property a' -<o o~ we conclude: 3rj G R, t G T 5(ir 3 -) =^o 
<t', HM(ir,) = HM(p'). So HM(/ij) ^ HM(p'), i = l,2, because in the case of equality rj would be 
signature-safe reductor for hi. It is possible only if HM(/ij) are canceled while subtraction with 
fc-coefficient, what means that HM(/ii) = HM(/i2). S-representation of hi is constructed by adding 
Kpo\y(h 2 ) to S-representation of (a',p r ). 
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Theorem 18. Every iteration of the algorithm after step \4a] the following invariant holds: for 
Ver eT,<T^ signatures of elements of B, exists r„ £ R, t a £ T : S{t a r a ) = a such that t a r a can't 
be signature- safe reduced by R. 

The set R a = {r £ R \ S(r)\o~} is not empty, because contains the element ro added during the 
first algorithm iteration with S(ro) = 1. Let r a be <H-minimal element of the set; take t a — s ° r ^ . 
Suppose that t a r„ can be signature-safe reduced by some r\ £ R. This gives that r% >h r a and 
both sides of inequality are non-zero. It means that during the iteration which inserts in R the 
last of {r a ,ri} the set B was extended by labeled polynomial t'r a , where t' — lcm ( h ^(^)>hm(tv)) 
and t'\t a . So we have iS^'r^) S(t a r a ) = a => S(t'r a ) =<; a -< signatures of elements of B. This 
signatures inequality implies that t'r a can't be element of B during the current iteration and was 
removed at the step|4d]of some previous iteration, so 3r2 £ Rr% <h t'r a ,S(r2)\S(t'r cr ). This is 
impossible, because the existence of r-i <h f a ,r2 £ Ra contradicts <H-minimality of r a . 

Theorem 19. Every iteration of the algorithm after step \4d\ the following invariant holds: V/i £ 
H,S(h) -< signatures of elements of B has S -representation over R. 

Suppose that invariant breaks during some algorithm iteration and take the -<o-minimal a 
that has non-empty corresponding set V a = {h £ H\h breaks invariant, S(h) = a}. Then R is 
S CT -basis. Mg £ I (0,g) has S-representation over {(0, gi), (0, g m )} C R, so a >~ 0. Select v a 
- one of the V a elements with -< q -minimal HM. It can't be signature-safe reduced by R because 
the reduction result v\ would be element oW a with HM(ui) -<q HM(w cr ). Take w a = t a r a from 
the invariant of theorem [TS] and apply lemma [T7] to v a , w a and R. The lemma says that v a has 
S-representation over RU All entries of w„ in the representation can be replaced by t a r a to 

acquire S-representation oiv a over R only. It's existence leads to contradiction. 

Lemma 20. If R is S-basis, then {poly(r) | r £ R} is a Grobner basis of ideal I. 

For \/p £ I we can take some h= (cr,p) £ H and apply lemma [T51 to it. 



Theorem 21. SimpleSignatureGroebner({gi, . . . , g m } , /) returns Grobner basis 



At the moment of algorithm termination B — so by theorem R is S-basis. 



Comparison with other algorithms 

The presented algorithm belongs to the family of the Grobner basis algorithms using signatures 
and being at some degree a modification of F5 algorithm from [3]. One of the main modifica- 
tion directions of F5 is simplifying and clarifying the connected theory usually bound with some 
thoughts about extending the area of inputs the algorithm can be applied to. Investigations in this 
direction can be found in [SJ|TTJ[TU]. The other direction is improving efficiency of computations 
by introducing criteria to detect and don't perform some unnecessary computations. It is studied 
in [5J [51 [3] and allows to perform computations in a way that reduces to the end only polynomials 
that are either new S-basis entries or corresponds to extending monoideal of signatures known to 
be zero polynomial signatures - called syzygy signatures. Generalization simultaneously using all 
criteria described in algorithms TRB-MJ and SB [9] achieve even more efficiency because all 
discardings are performed before any computational heavy operations like polynomial reduction or 
computation highest monomial of S-polynomial - so the non-trivial computations are never become 
unnecessary because their results are never discarded. 

All mentioned algorithms including original F5 use discarding criteria of two types: syzygy- 
based criteria and rewrite-like criteria, with separate proof of correctness for each type. The other 
common idea used in the algorithms are S-polynomials: even the algorithms that does not deal 
with S-polynomials directly make heavy usage of them in the correctness proof. 

This paper describes an algorithm computing minimal S-basis and discarding computations with 
efficiency identical to TRB-MJ, but using the only one discarding criteria in step !4dl which is based 



on <fj-ordering of R. The routine ReduceCheckingSignatures(cr, p, R) can use different reductor 



selection strategies and their efficiency is open question. The method proposed in this paper is 
based on the same ordering of R and is identical for homogeneous case with the methods used in 
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original F5. The correctness proof is given without the use of S-polynomials and is formulated in 
a way that allows to apply the clear algebraic interpretation of signature-based algorithms from 
[TT] to the presented algorithm. 

Algorithm simplification lead to simplification of programming its implementation and debug- 
ging. It is achieved by the smaller number of objects involved in computation and use of the 
same order for discarding criteria and reductor selection. Simplicity of implementation and the 
absence of complex data structures allows quick algorithm integration with any computer algebra 
system that can work with polynomials. The author's implementation linked below was written 
from scratch in a 8 hours what is a lot smaller than the time author spent implementing other 
algorithms with a similar tools. 

The algorithm was implemented in CH — h using low-level functions from computer algebra sys- 
tem Singular 3-1-4 and open source codes of C. Eder (one of the authors of [3]) for implementing 
F5-like algorithms in this system. The source is contained in "ssg" function in a file available at 
https://github.com/galkinvv/Singular-f5-like/blob/ssg/kernel/kstd2.cc 

Comparison of SimpleSignatureGroebner implementation with other Grobner basis algorithms 
implemented by C. Eder gives practical checks for the following theoretical facts: 

— algorithm SimpleSignatureGroebner correctly computes Grobner basis; 

— the number of polynomials in the result set is not greater than the number of polynomials in 
a result of other incremental algorithms that compute S-basis; 

— the execution time is not greater than execution time of other signature-based incremental 
algorithms. 
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